Efficient entanglement criteria beyond Gaussian limits using Gaussian measurements 
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We present a formalism to derive entanglement criteria beyond the Gaussian regime that can be 
readily tested by only homodyne detection. The measured observable is the Einstein-Podolsky- Rosen 
(EPR) correlation. Its arbitrary functional form enables us to detect non-Gaussian entanglement 
even when an entanglement test based on second-order moments fails. We illustrate the power of our 
experimentally friendly criteria for a broad class of non-Gaussian states under realistic conditions. 
We also show rigorously that quantum teleportation for continuous variables employs a specific 
functional form of EPR correlation. 
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Quantum entanglement plays a key role in making 
quantum-mechanical predictions distinct from their clas- 
sical counterparts. Its characterization and detection 
have been a topic of crucial importance that has attracted 
a great deal of effort, both theoretically and experimen- 
tally. Despite remarkable progress, there still exists a 
pressing demand for further developments particularly 
for continuous variables (CVs). In the CV regime, Gaus- 
sian entangled states have been the subject of study, 
but now considerable attention has been directed to non- 
Gaussian states. This may be attributed to the fact that 
non-Gaussian entangled states not only provide practi- 
cal advantages over their Gaussian counterparts, e.g. in 
quantum teleportation [l[ and dense coding 0, but also 
turn out to be essential, e.g. for universal quantum com- 
putation Q and nonlocality test 

The quantum state of a CV system is fully described 
in phase space by representing two quadrature variables 
X (position) and P (momentum). To address the cor- 
relation of two CV systems, the EPR operators v! = 
\g\XA — and v' — \g\PA + ^Pb can be denned, 

as Einstein, Podolsky, and Rosen (EPR) originally envi- 
sioned [5] . Here g is an arbitrary real number and A and 
B label each subsystem. When the sum of the variances, 
E[ = (A 2 u') + (A 2 {)'), is smaller than \ [g 1 + pj, the 

state of the total system is inseparable Q and such states 
are said to be EPR-correlated. In fact, it has been shown 
that the EPR correlation, after local unitary operations, 
is a necessary and sufficient condition for two-mode Gaus- 
sian entanglement n this Letter, we study a simple 
case, where we use {u,v} instead of {u',v'} assuming 
g = 1). Using the uncertainty principle as a requirement 
of a physical system in conjunction with the partial trans- 
position [7], Simon found a necessary and sufficient con- 
dition for Gaussian entangled states (Simon criterion) Q; 
this criterion is also concerned with the second moments 
of the quadrature variables |9j. These criteria are theo- 
retically easy to calculate and readily tested experimen- 
tally since the quadrature variables can be measured us- 



ing highly- efficient homodyne detectors. 

Beyond the Gaussian regime, some entanglement crite- 
ria have also been proposed I3"12|- Shchukin and Vogel 
particularly derived a hierarchy of entanglement condi- 



tions using higher-order moments [13l . Il4j . For finite- 
dimensional systems, a remarkable criterion based on co- 
variance matrices that unifies other criteria as its corol- 
laries has been developed (IBj . However, these tools 
are not necessarily practical for infinite-dimensional CV 
states because they may be nontrivial to_ assess theoret- 

Some crite- 



18] 



ically and implement experimentally 
ria are at least theoretically easier to assess but they do 
not seem useful for a large class of non-Gaussian states. 
Therefore, new entanglement criteria beyond the Gaus- 
sian regime must be developed, desirably in an experi- 
mentally friendly form, which is the objective of the cur- 
rent work. 

Here we present a formalism employing the usual EPR 



operator Oepr = u 



in an arbitrary functional form, 



which is readily testable with highly efficient homodyn- 
ers and is theoretically convenient to calculate. These 
functional EPR criteria turn out to detect CV entangle- 
ment well beyond the Simon criterion, thereby proving 
genuine non-Gaussian entanglement. Our criterion can 
be efficiently tested under realistic conditions including 
environmental interactions and detector inefficiency. We 
illustrate the power of our formalism by considering a 
broad class of non-Gaussian states and show that it can 
be a crucial tool in addressing CV problems of fundamen- 
tal importance, e.g. the robustness of Gaussian versus 
non-Gaussian entangled states against decoherence. 

In order to present the new entanglement criterion, we 
briefly revisit the so-called Braunstein-Kimble (BK) pro- 
tocol which is the standard for CV teleportation [19| . Un- 
like the original description, we study the protocol in the 
Heisenberg picture in order to gain a crucial insight. Let 
ain be an unknown input which is to be teleported, and 
a a and as the entangled quantum channel sent to Alice 
and Bob, respectively. Each field can be decomposed into 
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real and imaginary parts (two quadrature amplitudes), 
a, = Xi + iPi (i — in, A,B). In the BK protocol, Alice 
first superposes her mode with the input at a 50:50 beam 
splitter to obtain two output modes, a± = {a- la ±a,A) / 
and measures the X (P) quadrature of mode ct_ (&+)• 
The measurement outcomes {A_,P + } are transmitted 
to Bob who performs the displacement D(5) = e Sa ~ 5 a 
with 5 = \/2{X- +iP+) [Si- Although A_ and P+ are 
classical numbers, if formally represented by quantum 
operators [2l| . Bob's output becomes 



flout — o-b 



V2( 



X- + iP 4 



a\ + a B - (1) 



We rigorously justify Eq. (fTJ) by showing that it leads to 
exactly the same mapping in the Schrodinger picture [22 1 
(Supplemental Material). 

Teleportation fidelity — As a figure of merit to assess 
the performance of the BK protocol, the output fidelity 
F averaged over all input coherent states \jS) — D((3)\0) 
(/3: amplitude) can be used. Noting that the mean am- 
plitudes of input and output fields always agree in the 
BK protocol [20j, we realize that it suffices to look into 
the fidelity for the vacuum state input |0). Let Ut be 
the global unitary operator that gives the output state 
Pout = TrA,B{^t|0)(0|i n (8)pAB6 r /} after teleportation via 
a quantum channel pab (Supplemental Material). The 
fidelity is given by F = (0\p ont \0) = Tr{p out |0) (0|} = 
Tr{/9 out : e -a ° uta ° ut :} with the normal-ordered expan- 
sion of the vacuum |0)(0| =: e~ a+a : [23]. Using Eq. Q 



and the vanishing contribution of the input vacuum state 
under normal-ordering, the output fidelity is reduced to 
an average over the quantum channel /?ab , 



F = 



(-iy 



in 



out "out/ 



-A-'u-A'ti 



PAB 



(2) 



using the identity (a B — cla)(o.b — a a) — A 2 u + A 2 v 
[iol ]. We see now that all orders of EPR correlations 
E m = ((A 2 u + A 2 v) m ) contribute to the fidelity nontriv- 
ially, and this can be used to address various issues on 
quantum teleportation using non- Gaussian channels. As 
an example, due to an inequality e" a > 1 — a 2 , we de- 
duce F > I — Ei. Therefore, the condition Ex < i alone, 
using Gaussian or non-Gaussian channels, guarantees the 
fidelity over the classical limit | [24| . 

For a Gaussian channel, higher-order correlations are 
constrained by the lowest one £1, and E\ < 1 becomes 
a necessary condition for faithful teleportation. Surpris- 
ingly, using the formulation above, we find that for a non- 
Gaussian channel, it is possible to have near-unit fidelity 
even without the usual EPR correlation, E\ > 1. There- 
fore, the EPR correlation is not a necessary condition for 
a successful teleportation in the non-Gaussian regime. 
Remarkably, beyond the teleportation context, this sug- 
gests a possibility of detecting genuine non- Gaussian en- 
tanglement, for which the Simon criterion fails. Note 



that (i) the fidelity above 1/2 proves quantum entangle- 
ment [24[ and (ii) the fidelity can be predicted also by 
directly measuring E\j of the quantum channel over N 

lj . (See also 



runs, as F 



Implementation below and Fig. 1 (a).) 

Entanglement criteria — We present a formalism to 
construct an arbitrary functional form of EPR correla- 
tion as an entanglement criterion. Given a general func- 
tion we consider an ensemble average of the Hcrmitian 
operator P(Oepr), with Oepr = u 2 + v 2 - As detailed in 
Supplemental Material, for separable states, the value of 
(^(Oepr)) is bounded, 



•Fmin < (-F(Oepr)) < Tn 



(3) 



Here the upper (lower) bound for separability is given by 
-7 r max(min) = max(min)„{0„} over 71 = 0, 1, ... , where 

O n =4 dxdyT(2x 2 )e-?y 2 L n (y 2 )J (2xy)xy (4) 



(L n : Laguerre polynomial, Jo: Bessel function). Eq. (fj| 
is readily integrated for a regular function J- . If the in- 
equality (TTfl) is violated for any function J 7 , the state is 
entangled. 

Our method, when applied to the lowest-order EPR 
correlation E[, recovers Duan et aZ.'s criterion. For 
higher-order criteria, using T(z) — z m , we obtain sepa- 
rable conditions E m = ((A 2 u + A 2 v) m ) > ml. Below we 
adopt a general T in the Taylor-series form and demon- 
strate its usefulness in detecting CV entangled states 
broadly. 

Implementation — Separable states must satisfy the in- 
equality (|17[) also under any local unitary transforma- 
tions. This reasoning leads to a comprehensive set of 
entanglement criteria based on the EPR operator using 
generalized orthogonal quadratures. That is, the quadra- 
ture amplitudes under rotation, X[ — cos <piXi — sin^Pj 
and P[ = sincjiiXi + cos <piPi (i — A, B) can be employed 
and the violation of p7|) for any single pair {4>a, 4>b} 
is sufficient to verify entanglement. Given an entangled 
state, the violation usually occurs over a range of values 
cf>A and 4>b for a fixed J 7 , making an experimental test 
robust. 

In experiments, the generalized Oepr can be measured 
by first phase-shifting each mode with (pi (i = A, B) 
and injecting them to a 50:50 beam-splitter (Fig. 1 
(a)). One then measures the fixed quadrature X\ and 
Pi at each output, respectively, using homodyne detec- 
tors. [X' A - X' B = V2Xx and P' A + P' B = %/2P 2 ). We 
emphasize that neither a full state tomography nor the 
reconstruction of a probability distribution is needed for 
this test of entanglement: The EPR value Oepr.j ob- 
tained at each run (j = 1,...,N) is plugged into a 
test function T to construct the average over N runs as 
(P(Oepr)) = a? Y,j=i -^(Oeprj)- It is then compared 
with the bounds in Eq. (f?|) to verify entanglement. As 
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FIG. 1. (a) Experimental scheme to detect generic CV entan- 
glement by only two homodyne detectors (HDs). 4>a,b- phase 
shift, (b) (J-(Oepr)) — fma as a function of v (coherent am- 
plitude) and n (overall efficiency including transmission loss 
and detector efficiency) for a dephased cat state (p=0.3), with 
T(z) = e~ Cz (1+Dz). At each point {u, n}, the positive value, 
which verifies entanglement, was optimized over the parame- 
ters C G (0, fO] and D G [-80, 80]. (c) {T(6epr)} - J™ as 
a function of v (p = 0.5), with the error bar 8 e = Aj^/y/N in 
Eq. ©. (d) {LF{6 E PK))-^ mex )/6 e forastatec |00) + ci|ff) 
with T = e~ Cz (circle) and T = e~ Cz {l + Dz) (triangle). In 
(c) and (d), n = 0.7 and nth = 0.07 (electronic noise) with 
iV = 10 5 . 



one obtains a finite number N of data, it is of practical 
importance to see whether such a finite number of data 
suffices to beat the separability bounds properly. This is- 
sue can be addressed by looking into the statistical error 
S e of the mean, S e ~ Ajf / y/~N, where 



(^ 2 (Oepr)) - (^"(Oepr)} 1 



(5) 



We include this analysis below to demonstrate the power 
of our criteria under realistic conditions. 

Trial functions — For an entanglement test, we employ 
a function J-(z) = e~ Cz (l + J2m=i D m z m ^ , an expo- 
nentially decreasing function multiplied by a finite poly- 
nomial. This form can be regarded as the most gen- 
eral in view of the Taylor-series expansion of a regular 
function (M — > 00). Our criterion becomes more power- 
ful by adopting higher-order polynomials; Given an en- 
tangled state, our task is to find the set of parameters 
{C, ,D M } for which the inequality (fTT)) is vio- 

lated. If it occurs at a level M, so does it at higher levels 
M' > M. The same logic also applies to the degree of 
violation measured by, e.g., ((J"(Oepr)) - ^max)/4 op- 
timized over the parameters {C, D\, ■ ■ ■ ,Dm}- That is, 
the larger the M value, the higher the confidence level 
of the test for a given number N of homodyne mea- 
surements [Fig. I (d)]. In this Letter we illustrate our 



method by confining the polynomial to the first order, 
F(z) = e~ Cz (l + Dz). Remarkably, a broad class of 
non-Gaussian entangled states of current theoretical and 
experimental relevance can be detected with it. 
Using F{z) = e~ Cz (l + Dz), Eq. flU gives 



(1 - C)"- 1 
(1 + C) 



n+2 



[l-C 2 + D(l-C + 2n)] , (6) 



whose maximum is relevant for entanglement detection. 

Case (i) D = 0: Max {O n } is y^c when n = 
0. Therefore, if there exists any C > such that 
(e~ c ° EPR ) > jrjj, the state is entangled. For example, 
with C = 1, the separability bound becomes 1/2, which 
is known as the classical limit for teleportation [24J. 

Case (ii) D ^ 0: In this case, the number n yielding 
maximum O n in Eq. ([5]) depends on C and D. It should 
thus be obtained case by case, which is still straightfor- 
ward to calculate. For a given state, the violation of (fTT)) 
usually occurs in a broad range of C and D, which makes 
our scheme robust against noise. Below we consider pos- 
sible experimental imperfections like channel loss or ther- 
mal noise during signal transmission and detector ineffi- 
ciency or electronic noise in homodyne detection. These 
errors can be modeled by 



(i = A,B), 



(7) 



where a^d is the mode actually detected, di i0 the original 
mode, and i>i the vacuum or thermal mode characteriz- 
ing noise. The results below are all obtained analytically. 
Our method is also compared with other CV criteria di- 
rectly testable by homodyne detection such as the Simon 
criterion Q and the entropic criterion 16. 17|. 

Examples — Let us first consider a class of dephased 
cat-states [2^|, p cat ~ \v, v){v, v\ + | — v, — v, — v\ — 
p (\v, v)(— v, — v\ + I — v, — v)(y, v\), where \v) is a coher- 
ent state of amplitude v and 1 — p > represents the de- 
gree of dephasing. They are entangled for all v > and 
p > 0, which is not detected by the Simon criterion. The 
entropic criterion detects entanglement only for a large 
v and p, despite the improvement with the generalized 
Renyi entropy 3, Tt| • In contrast, our criterion detects 
entanglement for any v > and p > 0, with <pA = <t>B = 
in Fig. 1 (a). Considering further the overall efficiency 
rj characterizing channel loss and detector efficiency, we 
plot the value of (.F(Oepr)) — ^max in Fig. 1 (b), where 
each point is optimized over the parameters C £ (0, 10] 
and D £ [—80, 80] for a highly dephased state (p = 0.3). 
The positive values in all ranges imply the detection of 
non-Gaussian entanglement with the degree of violation 
increased with v and 77. Moreover, let us address a realis- 
tic situation where a total number N of data is obtained 
from homodyne detection with detector efficiency rj and 
electronic noise, which can be described by mixing signal 
with a thermal state (nth- thermal excitation number) 
[2^]. In Fig. 1 (c), the case of dephased cat (p = 0.5) 
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is plotted as a function of v for rj = 0.7 and n t h = 0.07 
(electronic noise), with the error bar. The positive values 
above zero clearly verify entanglement. 

Secondly, let us consider a state \^)b — co|00) +ci|ll). 
We have checked that the Simon criterion and our 
method detect entanglement for any cq and c\ values, 
whereas the entropic method .17] does so for only |co| > 
0.503. In Fig. 1 (d), we plot ((-F(O epr )) - F max )/5 e for 
the state | v E')b, which is improved by changing the test 
function from T — e~ Cz to T = e~ Cz (l + Dz). Fur- 
thermore, suppose that the state undergoes decoherence 
through a thermal reservoir. The reservoir interaction is 
described by Eq. (J7J). As rj gets larger the system be- 
comes more decohered. As shown in Fig. 2 (a), entangle- 
ment detection is significantly improved over the Simon 
method (square) , particularly by modifying the test func- 
tion from T(z) — e~ Cz (circle) to J-(z) = e~ Cz (l + Dz) 
(triangle). Similar results are also obtained for higher- 
dimensional entangled states. 

Finally, the classes of entangled states generated 
by practically available non-Gaussian operations, i.e., 
photon-subtraction (2|| and addition 27 1 are of current 
interest 



28|. 



When these operations are applied to a 
two-mode squeezed state |TMSS) = e s ( a ^-<^ aB )|00), 
the Simon criterion fails to verify the entanglement of the 
photon-added state a} A a} B |TMSS) for s < 0.378, whereas 
our method detects entanglement for the whole range 
of s > 0. Furthermore, we address a question of fun- 
damental importance raised recently by Allegra et al, 
i.e., whether Gaussian entanglement is more robust than 
non-Gaussian entanglement against decoherence due to 
the thermal reservoir j^]. We show that non-Gaussian 
entanglement can be more robust than Gaussian en- 
tanglement in the parameter regime they studied [30j |. 
The conclusion of [29( was drawn by using several previ- 
ously known criteria for verifying non-Gaussian entangle- 
ment, which however did not outperform the Simon crite- 
rion. In contrast, our method analytically shows that the 
photon-subtracted states (triangle: our criterion, square: 
Simon criterion) survive longer under thermal noise than 
their Gaussian counterparts of the same initial energy 
(circle) [Fig. 2 (b)]. Similar results are also obtained for 
the photon added states. 

Remarks — We have developed a method to derive en- 
tanglement criteria beyond the Gaussian limit that can 
be efficiently tested by homodyne detection. We illus- 
trated our method by adopting two trial functions which 
broadly detect non-Gaussian entanglements of funda- 
mental and practical interests. We also showed that 
our method is robust against experimental imperfections. 
The trial functions can be made more powerful by incor- 
porating higher-order polynomials, which do not however 
incur any extra experimental efforts. On another side, 
the question on how much useful (robust) non-Gaussian 
entangled states are under practical conditions is of im- 
portance for quantum information processing. Our study 
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FIG. 2. Entanglement is detected till the decoherence 
time log7? -1 for the initial state (a) co|00) + ci|ll) and 
(b) &a&b | TMSS) interacting with a thermal reservoir of (a) 
THh — 0.5 and (b) 0.05, respectively. Triangles: our method 
with J-(z) — e~ Cz (l + Dz), squares: Simon criterion, and cir- 



cles: (a) our method with T{% 



(b) Simon criterion 



for |TMSS) having the same initial energy as aAas|TMSS). 



for quantum teleportation, clarified here as employing 
the EPR correlation in a specific form, shows that non- 
Gaussian channels can be more robust against decoher- 
ence than Gaussian ones. 
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SUPPLEMENTAL MATERIAL 

Part I: We first show that Eq. (1) is a legitimate trans- 
formation describing the CV quantum teleportation that 
leads to exactly the same mapping in the Schrodinger pic- 
ture [22]. In Ref. [21], Eq. (1) was used for the study of 
teleporting temporal correlation functions by only noting 

that the operator 5 = \/2 (^X- + iP+^j = ai n — a A acts 

like a classical number due to [8, 8^] = 0. In contrast, we 
rigorously justify Eq. (1) as follows. 

Note that a out = y/2 ain +f- B — a A can ^ c understood 
as superposing two modes a; n and as at a 50:50 beam 
splitter followed by the nondegenerate parametric am- 
plification (gain G = 2) with the idler mode a a. This 
suggests the conjugate evolutions for a a and a B , 

a A =-\/laA 7= — , a B = 7^ — • (°) 



With the Heisenberg transformations in Eqs. (1) and 
(8) taken as = U}a,iU t {i — in, A, B), where U t is a 



global unitary operator corresponding to CV teleporta- 
tion, one can obtain the characteristic function of the out- 
put, C out (A) = tr{£>(A)/w}, where p out = tr A , B {U t p in ® 
PabU}}, as 

C out (A) = an(A)C AB (A*,A), (9) 

which is in agreement with [22]. 

Part II: We next prove that there exist quantum chan- 
nels that do not have the usual EPR correlation, i.e. 
Ei > 1, but can tclcport any input state with fidelity 
arbitrarily close to unity. For simplicity, let us consider a 
class of entangled channels produced by injecting identi- 
cal fields (with only a phase shift to the second) of modes 
a\ and di to a 50:50 beam splitter. The output modes 
a a = (ai + fl2)/v / 2 and clb = (—Si + 0-2) /V% give the 
fidelity in Eq. (2) as F = (e^ 2 ^ ) (e^ 2 ). For high fi- 
delity, we only need to maximize fi = (e~ 2Xl ) for a 
single mode; = (e~ 2P2 ) can be similarly maximized 
by applying | phase-shift (X — > P). To ensure the no- 
EPR-correlation, E\ > 1, for this channel, we impose the 
condition (X 2 ) = \ (no squeezing). Our task is now to 

make (e~ 2X i) close to unity with (X 2 ) = j. Obviously, 
the quadrature distribution p(X) must be sharply cen- 
tered around X = to obtain a large (e~ 2x2 ), but pre- 
sumably with a rather long tail to simultaneously ensure 
a finite variance, (X 2 ) = j. This leads us to envision 
a probability distribution that decreases with t = 2X 2 
more rapidly than the Gaussian function. We can con- 
struct a number of such distributions analytically, e.g. 

Pm(t) ~ e~ ctm (to > 1). After normalization, we obtain 
a class of probability distributions 

Pm (X) = ^\X\e-C-^)^, (10) 

and find that /1 = J dXp m (X)e~ 2x2 approaches unity 
as to increases. 

One can construct a corresponding non-Gaussian 
entangled state, p A B = UbsPi ® P2^b S (^bs: 
beam-splitter) with two single-mode states p\ = 
J dX Pm (X)\X)(X\ and ^ = / dPp m {P)\P}(P\, where 
\X) and \P) arc X- and P-quadrature cigenstates with 
the probability distribution p m in Eq. (10). The uncer- 
tainty principle yields A 2 P > 16 ^i x = \ for the state 
p\. No squeezing is then assured for arbitrary quadra- 
tures X cos 9 + P sin 9 with a little algebra. Further- 
more, it can be shown that the state pab asymptoti- 
cally gives the unit fidelity for all input states; Its char- 
acteristic function Cab(A*,A) = C p {\ x )C p {\ y ), where 
Cp(fc) = / dXp m (X)e^ 2 ^ 2lkx , becomes very flat over 
the entire range of A with to increased, which ensures a 
high-fidelity via Eq. (9). As an example, in Fig. 3, we 
plot the Wigner function of the teleported state for the 
input Fock state |1). The overall features are successfully 
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FIG. 3. Wigner functions W(a x ,a y = 0) for the input Fock 
state |1) and the teleported state (indiscernible with fidelity 
F w 0.992) by the quantum channel in Eq. (10) for m = 50 
(a). In (b) a rather long tail is shown in the output (upper 
blue curve), which does not significantly affect fidelity. 



teleported except for a relatively long tail in the output, 
whose effect on fidelity is however negligible. 

Part III: Finally, we show how the separability con- 
dition in (3) can be derived employing an arbitrary func- 
tion T. It is well known that for a separable state, the 
density operator under partial transposition (PT) is still 
positive-semidefinite [7], thus satisfying all requirements 
imposed on a legitimate quantum state. As demon- 
strated below, this leads to an upper/lower bound of the 
value (J""(0epr)) that a separable state can take, where 
Oepr = u 2 + v 2 = (Xa - X B ) 2 + (P A + P B ) 2 . For each 
J-, one can find out such bounds as prescribed below, the 
violation of which is a sufficient condition for entangle- 
ment. 

First, note that when a a and &b are input modes to a 
50:50 beam splitter [Fig. 1 (a)], the operators u = Xa — 
Xb and v — Pa + Pb correspond to two commuting 
observables s/2Xi and V2P 2 , respectively, at the output 
modes. From the viewpoint of the output state, it is 
thus obvious that the Hermitian operator J-(Oepr), for 
any function J 7 , can be regarded as a proper symmetric- 
ordered operator. Therefore, one can generally express 
the quantum average (.F(Oepr)) by means of two-mode 
Wigner distribution Wab(<x, P) as 

(•F(Oepr)) 



d 2 ad 2 l3W AB (a,l3)T [(a x — j3 x ) 2 + (a y + /3 y ) 2 ] , 

(11) 

where the EPR operator Oepr = (Xa — X B ) 2 + 
(Pa + Pb) 2 is replaced by the corresponding classical 
amplitudes (x/y: real/imaginary parts). Under PT, 
the Wigner function changes as WAB(& x ,a y , f3 x , /3 y ) — > 
Wab(&xi OLy, /3 X , ~ Py) [8] and Eq. (fTTj) becomes 

M= (-F(Oepr))pt 

= ( d 2 ad 2 pW AB (a,p)F[{a x -P x ) 2 + (a y -p y f] . 



M = J d 2 a 1 d 2 a 2 W l2 (a u a 2 )J : (2\a x \ 2 ) (13) 

where Wi 2 (ai,a 2 ) = W A b ^f 2 ') represents a 

certain two-mode Wigner function. We now see that the 
extremal values of M. for all quantum states are the same 
as those of a single-mode PFi(ai) = J d 2 u 2 W\ 2 (ct\, a 2 ) 
as T now involves only one mode a\ in Eq. (|13p . From 
now on, we thus work with only a single-mode state p 
having a Wigner function W(a). Writing VF(a) as the 
Fourier transform of the characteristic function C(A) = 



Tr 
tain 



pD(X) , i.e., W(a) = ± J d 2 XC(X)e 



M 



d 2 aW(a)F (2\a\ 2 ) = Tr pO 



K , we ob- 



(14) 



where 



d 2 ad 2 \e aX "- c 



"D(X)T (2\a\ 



(15) 



Now the problem is reduced to that of obtaining extremal 
expectation values of the operator O in Eq. (fT4|) . Obvi- 
ously, the maximum (minimum) value of Tr pO is the 

maximum (minimum) eigenvalue of O. One can show 
by a direct calculation that the operator O does not 
have off-diagonal elements in the number-state basis, i.e. 
(m|0|n) = for m ^ n. Thus, the eigenvalues of O are 
the diagonal elements O n = (n\0\n) given by 

/>oo 

O n =4 dxdyT(2x 2 )e-i y2 L n (y 2 )J (2xy)xy (16) 



using (n\D(X)\n) — e~5l A l 2 L„(| A| 2 ). (L n : Laguerre 
polynomial, Jo: Bessel function) The ensemble average 
A4 = (O) in Eq. ([M)) can thus take a value only in the 
range [J" mi „, J" ma x] where -7 r mi n(ma X ) = min(max)„{0„} 
over n = 0, 1, • • • . Therefore, for any function J 7 , a sepa- 
rable state must satisfy 



(17) 



As an example, our method, applied to the lowest- 
order EPR correlation E[ = (A 2 u') + (A 2 v'), recovers 
Duan et aUs criterion. For this case, the argument of T 



inEq. (JEJ becomes (on-fii) 2 -> (\g\a.i 
Then, by redefining a single-mode a.\ = 



(i = x,y). 
f(\9\a- 



-0), we similarly obtain the bounds, Eq. (|16p. with only 



the replacement F(2x 2 ) — > J- 



^)x 2 
g 2 ' 



Identifying 



With the change of integration variables (redefining new 
modes) as ct\ = and a 2 = ^=r, the above expression as E[ = (A 2 u') 



(12) F(z) — z, one obtains the minimum O n — \(g 2j r^i) with 
n — 0. Therefore, we obtain the Duan et al.'s criterion 



(A 2 v')>\(g 2 + j,). 



